Fermi Coordinates of an Observer Moving in a Circle in Minkowski Space: 
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Coordinate transformations are derived from global Minkowski coordinates to the Fermi co- 
ordinates of an observer moving in a circle in Minkowski space-time. The metric for the Fermi 
coordinates is calculated directly from the tensor transformation rule. The behavior of ideal clocks 
is examined from the observer's reference frame using the Fermi coordinates. A complicated relation 
exists between Fermi coordinate time and proper time on stationary clocks (in the Fermi frame) 
and between proper time on satellite clocks that orbit the observer. An orbital Sagnac-like effect 
exists for portable clocks that orbit the Fermi coordinate origin. The coordinate speed of light is 
isotropic but varies with Fermi coordinate position and time. The magnitudes of these kinematic 
effects are computed for parameters that are relevant to the Global Positioning System (GPS) and 
are found to be small; however, for future high-accuracy time transfer systems, these effects may be 
of significant magnitude. 



I. INTRODUCTION 



> 

During the past several decades there have been significant advances in the accuracy of satellite-based navigation 
and time transfer systems [0. The most recent example of such a system is the Global Positioning System (GPS), 
which is run by the U. S. Department of Defense @-|§. The GPS is based on a network of 24 atomic clocks in space, 
from which a user can receive signals and compute his position and time. The order of magnitude of the accuracy 
of the GPS currently is ten meters and twenty nanoseconds ||. Future navigation and time transfer systems are 
"S^ • being contemplated that will have several orders of magnitude better position accuracy and sub-nanosecond time 
Q-i accuracy j6j-||]. 

Over satellite to ground distances, precise measurements must be interpreted within the framework of a relativistic 
theory The design of high-accuracy navigation and time transfer systems depends on detailed modeling of 
the measurement process. The comparison of experiment with theory in general relativity is more subtle than 
in special relativity or other areas of physics for two reasons. First, due to the possibility of making coordinate 
gauge transformations, often it is not clear what is a measurable quantity and what is an artifact of the coordinate 



system [ [10|J11|| . Second, all measurements must in-principle be interpreted in the comoving frame of the experimental 
apparatus [12]. In the case of the GPS, even theoretical relativistic treatments do not interpret the measured observable 
in the comoving frame of the apparatus, see for example Ref . fl3| . 

The principle reason for investigating in detail relativistic effects is to improve the current accuracy of GPS and 
to create future time transfer and navigation systems that have several orders of magnitude better accuracy. At the 
present time, it is well-known that small anomalies exist in position and time computed from GPS data. The origin of 
these anomalies is not understood. In particular, GPS time transfer data from the U. S. Naval Observatory indicates 
that GPS time is periodic with respect to the Master Clock, which is the most accurate source of official time for the 
U. S. Department of Defense [pT[-[l^]. Furthermore, other anomalies have been found in Air Force monitor station data 
that are not understood at present ]l^] . The work below is a first step in revisiting from a first principles approach 
some of the physics that has gone into the design of the GPS . 

A theoretical framework, known as the tetrad formalism, has been developed for relating measurements to relativistic 
theory |12 T^-|20|. However, due to a historical lack of accurate measurements over satellite to ground distances, in 



most cases it has not been necessary to apply the tetrad formalism in the area of satellite-based systems. This 
situation is rapidly changing due to improved measurement technology P,p3|-p5[]. 

In the tetrad formalism, measurements are interpreted as projections of tensors on the tetrad basis vectors. Con- 
sequently, these projections are invariant quantities (scalars) under transformation of the space-time coordinates. 
However, these projections do depend on the world line of the observer and the choice of local Cartesian axes used by 
the observer |2^,|l^] . The need for the tetrad formalism to relate experiment to theory, and the problem of measurable 
quantities in general relativity, are well discussed by Syngc Q, Soffel [EM and Brumberg |22]. 
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The tetrad formalism was initially investigated for the case of inertial observers, which move on geodesies 
jUI ^0| , p6| -p^]. However, many users of geodetic satellite systems, are terrestrially based observers, or are based 
in non-inertial platforms, such as missiles. The theory for the case of non-inertial observers has been investigated 
by Synge |Q2], who considered the case of non-rotating observers moving along a time-like world line, and by others 
p9|-|33||, who considered accelerated, rotating observers. For arbitrary observer motions, the effects are indeed very 
complicated, and the general theory gives limited physical insight. For this reason, in this work, I work out in detail 
a simple model problem. I treat the particular case of an observer moving in a circle in Minkowski space. I use the 
tetrad formalism [fl2| to obtain the transformation from an inertial frame in Minkowski space, to the accelerated frame 
of reference of an observer moving in a circle. The observer Fermi- Walker transports his tetrad basis vectors. The 
case presented here corresponds to an observer kept in a circular trajectory by rockets, i. e. , he is not on a geodesic. 
The explicit results presented below form a particular case of the general results obtained by Synge jl2| , specialized to 
zero gravitational field and circular motion. The situation treated here also corresponds to a reference frame that has 
a constant orientation with respect to the distant stars |34|, such as an Earth-centered inertial reference frame used 
in GPS satellite computations. The results are given here as expansions in small velocity of the observer, compared 
with speed of light. However, the coordinate transformation equations are not expanded in a power series in time, 
and hence, they are valid for arbitrary long times. 

This work explores the apparent behavior of ideal clocks as observed from a Fermi coordinate frame of reference of 
an observer moving in a circle in flat space-time. The Fermi coordinate frame is the closest type of non-rotating frame 
to that which is used in experiments and in GPS. In particular, a Fermi coordinate frame is a non-rotating frame of 
reference, where Coriolis forces are absent, and light travels essentially along a straight line [L2[. The observer's Fermi 
frame in this work is analogous to the Earth-centered inertial (ECI) frame used in GPS. The ECI frame moves with 
the Earth but maintains its orientation with respect to the distant stars. 

This work neglects all gravitational field effects that are significant for the actual behavior of GPS clocks. As such, 
this work treats only the kinematic effects associated with the circular motion of the Fermi coordinate frame about 
a central point (analogous to the sun's position). Since gravitational fields are neglected in this work, the observer's 
tetrad is Fermi- Walker transported along a time-like world line (a circle in 3-d space). In the case of GPS, the origin 
of the ECI frame follows a geodesic path (the Earth's orbit). Therefore, the observer's Fermi frame in this work is an 
accelerated, non-rotating frame, while the ECI frame is in free fall. As remarked above, the effects on clocks are in 
general very complicated, even without the gravitational field. For this reason, I work out in detail the specific case 
of an observer in circular motion in Minkowski space-time, and I apply the results to parameters that are relevant to 
the GPS. 



II. FERMI- WALKER TRANSPORT DIFFERENTIAL EQUATIONS 

Consider a time-like world line C of an observer moving in a circle in 3-dimensional (3-d) space given by 

x°=x% + u (1) 

x = x Q + a cos (uiu/c) (2) 

x 2 = Xq + a sin (uju/c) (3) 

x 3 = x 3 (4) 

where Xq, xj, Xq, and Xq are constants, u is a parameter along the world line, and c is the speed of light in an inertial 
frame in vacuum. The observer travels in a circle in the x — y plane. See Figure |l|. Along his world line, the observer 
travels with a four velocity 

u =it; (5) 

and has an acceleration 

where Y % - k is the affine connection. From the normalization of the 4- velocity, u % Ui = — 1, I obtain the relation between 
the arc length, s = ct, where t is the proper time, and the parameter u, 

u = 7s (7) 
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where 7 = (1 — v 2 ) 1 I 2 . The parameter v is the dimensionless velocity, given by 

au> 



(8) 



I use the convention that Roman indices take the values i = 0, 1, 2, 3 and Greek indices take values a — 1,2, 3, see the 
Appendix for the conventions used in this work. 

In Minkowski space, the affine connection components all vanish, T l - k = 0. Hence, the explicit expressions for the 
4- velocity and acceleration are: 

i* s =7 U. sm( ), — cos( ),0 (9) 

V c c c c / 



w l = —a j 2 



(^ 2 ( ,cos(^ S ),sin(^ S ),0) (10) 



As the observer moves on the time-like world line C, he carries with him an ideal clock and three gyroscopes. I 
use the notation x l for Minkowski space coordinates. At some initial coordinate time x°, the observer is at point Pq 
at proper time r = 0, corresponding to s = 0. On his world line, the observer carries with him three orthonormal 
basis vectors, A^, where a — 1,2,3, labels the vectors and i = 0,1,2,3 labels the components of these vectors in 
Minkowski coordinates. These vectors form the basis for his measurements [Q. See Figure |^. The orientation of 
each basis vector is held fixed with respect to each of the gyroscope's axis of rotation |55|. The fourth basis vector is 
taken to be the observer 4- velocity, A/qn = v\ The four unit vectors A/*\ form our tetrad jjjj, an orthonormal set at 
Pq, which satisfy 

9ij A( a) X\ b) =$ab (11) 

where the Sij = 1 if i = j, and zero otherwise. 

At a later time s = cr, the observer is at a point P. The observer's orthonormal set of basis vectors are related to 
his tetrad basis at Pq, by Fermi- Walker transport. Fermi- Walker transport preserves the lengths and relative angles 
of the transported vectors. Furthermore, Fermi- Walker transport of basis vectors along the observer's world line 
corresponds to keeping the vectors as parallel as possible by keeping their orientation constant relative to the spin 
axes of the gyroscopes ]3^ |. For a given vector, with contravariant components p in Minkowski space, its components 
at P are related to its components at Po by the Fermi- Walker transport differential equations fl2|| 

SP 



where 



W i] = u l w l - w l u j (13) 



When we use Eq. (|l^) to transport a vector /* that is orthogonal to the 4- velocity, u 1 ji = 0, the second term in Eq. jl^ ) 
does not contribute. We refer to transport of such space-like basis vectors as Fermi transport, and — > = u l w l . 
For an arbitrary vector f l perpendicular to the 4- velocity, the explicit form of the Fermi transport differential equations 
for the world line in Eqs. (0)-(§) is 

df° 

= - 7 3 ^ cos ( 7 f 1 - 7 3 ^ sin ( 7 £) f 2 (14) 

■j- = 7 3 ^ 2 cos ( 7 sin ( 7 f 1 + 7 3 ^ 2 sin 2 ( 7 f (15) 
df 2 

At = "7 3 ^ 2 cos 2 ( 7 f 1 - 7V cos ( 7 sin ( 7 f (16) 
<K 

df 3 , , 

-j- = 17 
d£ 

where the components p are functions of ^ and the dimensionless proper time is given by 
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e = -- as) 

c 

The differential equations in Eqs. (p^)-(p7|) are nothing more than the equations that describe the Thomas precession 

of an electron's spin vector as it moves in a circular orbit around the nucleus The solution of Eq. (|l4|)-(|l7j) is 
given by |$5| 

f = -~tvA cos ( 7 2 £ + a) + (3 (19) 

f 1 = A cos (7£) cos (7 2 ^ + a) + A-f sin (7^) sin (j 2 ^ + a) (20) 

f = A sin (7^) cos (-f 2 £ + a) - Aj cos (7^) sin (j 2 £ + a) (21) 

f = 5 (22) 

where A, a, 5 and f3 are integration constants. 



III. DETERMINATION OF THE TETRAD 



At the point Pq, the tetrad basis vectors satisfy Eq. (|l l|) . These equations constitute 12 relations (since Eq. (11) is 
symmetric in a and (3) for 16 components A| a \ . However, the zeroth member of the tetrad is the 4- velocity, which is 
assumed known 

A{ 0) =u*= 7 (l,t;\t;V) (23) 
and is given by Eq. (g). Equation (0) can be rewritten as two equations 

S«4oA(0) = 0, a =1,2, 3 (24) 

9ii*U X U =Sab > a=1 > 2 > 3 (25) 

Substituting the explicit form of the Minkowski metric (see the Appendix) into Eq. (|2J), gives a condition on the 
zeroth components of a tetrad vector if the spatial components arc known (Greek indices take values a,K =1,2,3) 

A° Q) = S^vPXfa (26) 



Then use of Eq. (26) to eliminate the time components \9 from Eq. (|25|), leads to an orthogonality relation for only 



the spatial components of the tetrad vectors 

<W A (a) A 03) - 5^S K7 v fl v K X 1/ ia) X y {[3) = 5 a p (27) 
The general idea is to solve Eq. (|27|) for the spatial components of the tetrad at s — 0, and then to substitute the 



spatial components into Eq. (26) to obtain the time components of each tetrad vector. Thus, I obtain the tetrad 
vectors at the initial time corresponding to s = 0. For each tetrad vector, these components are then used as initial 
conditions in Eq. (||^ )to determine the tetrad components at point P at finite s. 

The exact solution of Eq. ( p7j ) can be obtained by noting that for v a = 0, the spatial components of the tetrad 
vectors are orthonormal, and the solution is given by 

*U = S£ (28) 

I assume that the tetrad unit vectors are approximately oriented parallel to the observer's local Cartesian x, y and z 
axes. For small v a << 1, the solution of Eq. ( p7| ) can be obtained as a triple power series in v a , by iteration, starting 
with Eq. (p8| ) as the first approximation. Having solved for the tetrad components as a power series, it is easy to 
guess the exact solution to be 

A? a) = 7<W (29) 

X U = *S+ 2 ^*«a«*w" (30) 

A l (0) =n i = 7(l,« 1 ! « 2 ^ 3 ) (31) 



4 



where the components of dimensionless velocity v a are given in Eq. (||). 

Equations (|29|)-(|3l|) give the components of the observer's tetrad basis vectors at point Pq. At point P, the arc 
length is s > 0, and Eqs. (p9|)-(31 ) give the initial components for each tetrad vector. These initial components are 
used in the general solution of a Fermi- Walker transported vector, given in Eq. (0). The tetrad components at point 
P at finite s are given by 



A 



(0) 



{7, -7^sin(7£), 7^cos(7^), 0} (32) 
A (i) = {"7^ sm(7 2 0,cos( 7 £) cos(7 2 +7 sin( 7 sin( 7 2 f), cos( 7 2 f) sin( 7 f) - 7 cos(7 £) sin( 7 2 £),0} (33) 
*(2) = {7^ cos(7 2 0,-7 cos( 7 2 Bin(70+oos(70 sin( 7 2 £),7 cos(7f) cos( 7 2 f) + sin( 7 £) sin( 7 2 £),0} (34) 



A^ 3) -{0,0,0,1} 



(35) 



Equations (|3^)-(|35|) are valid for all times s because they are exact solutions of Eqs. Jl4|)-(|l7|) and satisfy the 
orthogonality relations in Eq. (fi"i"|). 



IV. FERMI COORDINATES 

Fermi coordinates are defined by the geometric construction shown in Figure ^. Every event P' in space-time has 
coordinates x 1 ' in the global inertial coordinate system. According to the observer moving in circular motion along 
a time-like world line, the same event has the Fermi coordinates X^ a \ a = 0,1,2,3. The first Fermi coordinate, 
X(°> = s, is just the proper time (in units of length) associated with the event P'. The proper time for P 1 is defined 
as the value of arc length s in space-time at which a space-like geodesic from point P (which defines s) and passes 
through event P' . The definition of s requires the tangent vector of this space-like geodesic, /i*, be orthogonal to the 
observer 4- velocity at P, 

^Vi\ P = (36) 

For the simple case of Minkowski space, the geodesic is simply a straight line connecting the points P P' 

// = N(x u - x l ( 1S )) (37) 

where x*(u) and x %l are the coordinates of P and P', respectively, in the global inertial frame, and A is a normalization 
constant that makes fi 1 a unit vector. The orthogonality condition in Eq. ([56]) is 



9ijtf(s) A ( )( s ) = 



(38) 



and gives an implicit equation for s for a given P' . This orthogonality condition gives the first Fermi coordinate of 
the point P' 



X (0) 



and Eq. (|3q) gives s as an implicit equation: 



7 s = x 



(x L - Xq) sm(- ) - {x zl - Xq) cos(- ) 



(39) 



(40) 



In the limit of small speeds of the observer, v — > 0, Eq. 
for s by iteration, resulting in a power series in v 



c c 
lives s — x° ' — Xq . For v « 1, Eq. (fib]) can be solved 



X<°> = 8 = Ax° 
1 

~2^ 



Ax sm( ) - Ax cos( ) 

c c 



v Ax° v 2 



-a + 2 Ax 1 cos( ) + 2 Ax 2 sin 

c c 



Ax 2 cos(- 



) - Ax 1 sm( ) 



v^ + Oivf (41) 



where I use the notation 
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/^X — X 



(42) 



For future reference, I label the coordinates of P, P', and Pq in the global inertial frame by P' — (a; ' ,x 1 ' ,x 2 ' ,x 3 '), 
P = [x , x , x , x ), and Pq = (xq, x , x q , x ). 

The contravariant spatial Fermi coordinates, X^ a \ a = 1,2, 3, are defined to be jl2| 

*W = a ^ a) = 9ij a(s) J (s) r,™ % (s) (43) 

where I used the Minkowski definition of space-time distance, a, along the space-like geodesic between P and P' 



o*=g ij (x , '-x t ){r" 



(44) 



where gij is the Minkowski metric and rf- 1 ^ = fjuj) = Vij i s an invariant matrix, which is also equal to the Minkowski 
metric, given in the Appendix. The parameter s in Eq. ( J4^ ) is a function of P' coordinates x 1 ' , so s must be 
eliminated using Eq. (j4l|), so that are explicit functions of the global inertial coordinates x 1 ' . I have calculated 

this transformation correct to fourth order in v. However, since the expressions are complicated, I write them here 
correct only to third order in v. 



I' 1 ' = Ax 1 - a cos( ) + - 

c 4 



Ax 1 - Ax 1 cos( 



Ax 2 sm( ) 



1 

2a" 



Ax° 



. 2 . ,Ax°uj' 

-Ax + a sm( ) 

c 



^ + 0(i/) 4 



X {2) = Ax 2 - a sm( ) + - 

c 4 



2Ai u w 



Ax 2 + Ax 2 cos( ) - Ax 1 sin 



-Ax° 
2a 



. , .Ax°a;.' 

Ax — a cos( ) 

c 



*/ 3 + O(^) 4 



X^ = Ax 3 

The inverse transformation can be obtained from Eqs. (^5|)-(|47]) by iteration resulting in: 
x Q/ ~x =X^ + 

,=XM+acos(^^) + 



XV c 0S (^)-X« sin(^) 
c c 



„1 / 



,uXW, 1 



X (L > - X (L > cos( ) - X (Z > sin 

c c 



X {2) +a sin + - 

c 4 

X (3) 



X( 2 > + X< 2 ) cos(^^) - X« sin(^^) 



^ 2 + 0{vf 



v 2 + 0(vf 



(45) 

(46) 
(47) 

(48) 

(49) 

(50) 
(51) 



V. METRIC IN FERMI COORDINATES 

The line element in the Fermi coordinate system of the observer is written as 

ds 2 = -G {ij) dX {i) dX {3) (52) 

where Guj\ are the metric tensor components when the X^ are used as coordinates. A direct calculation of the 
metric tensor components from the tensor transformation rule 

dx^ dx^ 

G m = 9kl (53) 

using the transformation Eq. p5|)-(^7|) and the Minkowski metric for g^, gives 

G ( oo) = -(l + C) 2 (54) 



G 



and G( Q/3 ) = 5 aj3 where ( is given by 

C 



Ax 1 ,wAj°. 
1 cos( ) sm( 



Ax 2 



loAx 



o 1 



^ 2 + 0(u 4 ) 



a c a c 

Using the transformation Eqs. d45|)-([47|) to express £ as a function of the Fermi coordinates, I obtain: 

1 



c 



X« cos(-X(°)) + sin(-X(°)) v 2 + 0{v A ) 
c c 



(55) 



(56) 



The result in Eq. (|I] ) agrees with Synge's general theory |12| of Fermi coordinates specialized to flat space. Synge 
shows that 



C = X (0) W W = V(af3) W W = v ia0 ) w l Ap } 



(57) 



To third order in v, Eq. (56) agrees with the general theory given by Synge |l^]. The quantities X^ — X^ are 

the covariant Fermi coordinates and — w l are the components of the observer's 4-acceleration in the Fermi 
coordinate system. The 4-acceleration, w l , i = 0, 1, 2, 3, is related to the ordinary acceleration, a 13 , (3 = 1, 2, 3, by 



w 



7 2 d 2 x l 
~^ 2 ~d£ 2 



7 



■ a 



c at* c* 

Using this relation, ( can be written in terms of the global inertial coordinates of P and P' as 

7 2 



C= i T S a p (x a '-x a (s)) a? 



Cf 



(58) 



(59) 



For the circular motion treated here, from Eq. (pq), the three spatial components of the observer's acceleration in 
Fermi coordinates are given by 



w 



(/3) 



V 

a 



<-X(°)),sin(-X(°)),0 



(60) 



where we have dropped fourth-order terms in v. To third order in v, the observer's acceleration is just the classical 
centripetal acceleration for an observer moving in a circle. 

The quantity c 2 £ is a potential (analogous to a gravitational potential) that determines the rate of proper time 
with respect to coordinate time in Fermi coordinates. Proper time is kept by an ideal clock. Coordinate time depends 
on the definition of the reference frame and coordinates used within that frame. The quantity £, given in Eq. (p3S|), 
depends on acceleration , and the difference in coordinates, x al — x a (s), from the spatial origin of the Fermi 
coordinates. Away from the spatial origin of Fermi coordinates the acceleration produces an effective gravitational 
potential field, c 2 £, which leads to a change in the relation of coordinate time to proper time, through the metric 
components G(y). 



VI. APPARENT BEHAVIOR OF CLOCKS IN FERMI COORDINATES 



The significance of proper time is that it is the quantity that ideal clocks keep, that it may be related to time 
on real clocks, and that it is closely related to measurable quantities. The significance of coordinate time is that it 
enters into the theory. In order to analyze experiments with clocks, we must relate what is measured, i. e. proper 
time intervals between space-time events, to coordinate positions and times of these events. 

Consider an ideal clock whose world line is given in Fermi coordinates by X^ 1 ' — jW (X^). At coordinate times 

X^ \ and Xg \ the clock is at spatial positions X^\ and X^\ respectively. The elapsed proper time between these 
two events, At — As/c, is given by the integral along the world line of the clock 

1/2 



As = 
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A. Stationary Clock 



For the case of a stationary clock, located at constant Fermi coordinate position X^ a \ the relation between proper 
time r and coordinate time, X^\ is given by taking dX^ = in Eq. (|f). The rate of proper time with respect to 
coordinate time is then given by 



ds 



i + Cpf (0) ,^ (Q) ) 



(62) 



Equation (|62j) gives the well-known result that the rate of proper time depends on the location of the clock [j42| . For 
the simple case of circular motion of the observer with ( given by Eq. (|56[) , substitution in Eq. ( |62| ) and integration 
from XjP to Xg^ gives 



As = X 



(0) 



X 



(o) 



c 



{XW [sin (^»)- sin (^°>) 



X (2) 



(^-eos^f)]} 



(63) 



where I have dropped terms of 0{v 4 ). Equation ( |63] ) shows that there is a complicated relation between elapsed 
proper time on a clock located at Fermi coordinates X^\ X^ 2 \ and Fermi coordinate time, X(°\ For example, 
a clock on the X^ axis (with X^ = 0) has a periodic time with respect to Fermi coordinate time, given by 



As = Xg — (au/c)X^ sm(ujX B ' /c). The initial conditions in Eq. ( |63| ) determine complicated phase relations 
between proper time and coordinate time. 



B. Clock in Circular Motion 





= fecos 












= &sin | 


^ x (o) 








X (3) 


= 





Now consider an ideal clock that is a distance b from the origin of Fermi coordinates and the clock executes circular 
motion in Fermi coordinates in the X^ — X^ plane. I take the world line for this motion to be given by 

(64) 

(65) 
(66) 

The clock moves in the X^-X 2 ^ plane, which coincides with the plane defined by the circular motion of the origin of 
Fermi coordinates. See Figure |[ The situation considered here is a simplification of the problem of a GPS satellite 
(with a clock) orbiting the Earth, which is itself orbiting the sun. (In the case of GPS, the satellites are in orbits that 
are inclined at 55° with respect to the Earth's equator.) However, as mentioned previously, in the case of the Earth 
and GPS satellite, both travel on geodesies. Since we are neglecting the effects of gravity in the problem considered 
here, the Fermi observer and orbiting clock are both Fermi- Walker transported, e. g., with the aid of rockets. 

Substitution of the satellite clock's world line into Eq. (|6l| ) , leads to the following relation between elapsed proper 
time observed on the clock and Fermi coordinate time 



As 



1 



1 (bCl 

2 \~c~ 



(X 



(o) 

B 



Xf) 



ab U3 



A B 



■v „(0) 
A A 



(67) 



The rate of proper time with respect to coordinate time is 

2 



ds 



dX(°) 



1 nn 

2 It 



abcos X y ' 

c / V c 



(68) 



The proper time has a constant rate offset from coordinate time, represented by the first term in Eq. (|67j). In 
addition, the proper time has a periodic component with respect to coordinate time, given by the second term in 
Eq. @. 

To demonstrate the magnitude of these effects, I use parameters that correspond to Earth-Sun distance and GPS 
satellites. I compute the magnitudes of the following terms with the values shown in Table |: 



1 An 

2 \ c~ 



8.348 x 1CT 11 



(69) 
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1 ( b Jl ) x ^ = 3.597 x 10- 6 sec (70) 

2 \ c J Si 

Equation ( |69| ) gives the constant rate offset of the moving clock with respect to Fermi coordinate time, due to time 
dilation. In GPS, this term leads to a net slowing down of GPS clocks by approximately 7 /is per day |37|| , which 
results from multiplying the value given in Eq. ( |70"| ) for a single revolution by a factor of 2 (since the GPS satellites 
make approximately two revolutions per Earth day.) 

The second term in Eq. ( |67j ) is a periodicity of the proper time with respect to Fermi coordinate time with amplitude 

nh /,i 2 

= 12.05 x 10~ 9 sec (71) 



c Q — i 

which corresponds to an amplitude in the rate 



2 = 1.755 x 10~ 12 sec (72) 



The proper time given by Eq. (|67]) is periodic with respect to Fermi coordinate time with the difference frequency, 
f2 — oj, which is the difference frequency between satellite and observer rotations. This is a kinematic effect due to 
the type of coordinate system used (Fermi coordinates) and the fact that the observer (Fermi coordinate frame) is 
moving along an arc of a circle, and therefore experiences an acceleration. Note that the periodic term vanishes when 
the observer angular velocity w = 0. If the observer were moving a straight line and hence have zero acceleration, 
w l = 0, the periodic effect would also be absent, as can be seen from Eq. (p7|). 



C. Sagnac-Like Effect for Portable Clocks 



Sagnac |3q-|40|] demonstrated that there exists a phase shift between two counter propagating light beams on a 
rotating platform and that this phase shift depends on the angular frequency of rotation of the platform. In the 
Fermi coordinates of the observer moving in a circle, there exists a Sagnac-like effect for two portable clocks that are 
orbiting in opposite directions. 

Consider a clock that moves eastward along a circle of radius b in the equatorial plane of the Fermi coordinates 
with world line given by Eq. (|64|)-(|66|). See Figure ||. The speed of the clock in Fermi coordinates is bfl/c. Assuming 
that the clock begins its journey at point A at time — and ends it at — 27rc/Sl, which corresponds to one 



revolution, the proper time that elapses on the clock is given by Eq. (67) as 
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The first term in Eq. ([73]) is the time dilation due to the speed bQ/c of the clock. The second term depends on the 
path traversed. 



Next, consider an identical clock that starts at point A at time X 



(0) 



but moves westward. The world line of 



this clock is then given by Eq. fl64|)-(| 
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|fi|. The proper time on this clock is given by 
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When the east-moving and west-moving clocks have each made one revolution and they have returned to their starting 
positions at point A, the difference in their times is given by Ar + — At_ = (As+ — As_)/c. For the case applicable 
to GPS, 6 << a and uj << |S1|, the difference on the two clocks is given by 
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(75) 
(76) 



where in the last line I have dropped terms of 0(w 2 /Sl 2 ). This Sagnac-type effect is due to the non-inertial nature of 
the Fermi coordinate system. The origin of this effect can be understood from the point of view of the global inertial 
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coordinate system. If the Fermi coordinate origin moved along a straight line in Minkowski space (rather than a 
circle) then there would be no acceleration and C = 0. In this case, the east-moving and west-moving clocks would 
have zero time difference after one revolution. However, since the origin of Fermi coordinates moves along a circle, the 
acceleration picks out a direction in the Fermi coordinate space-time and the symmetry between the two clock world 
lines is broken. Consequently, the elapsed proper time is different on the two clocks when they are compared after 
one revolution. Note that the time difference on the two clocks increases with increasing angular velocity of rotation 
u>, analogous to the optical Sagnac effect [B8H40L This increase is due to the fact that the breaking in symmetry is 
due to the magnitude of the acceleration, w^^see Eq. (|57|). Using the values of the parameters appropriate to GPS 
satellites, given in Table 1, I find the magnitude of the time differernce on the two clocks after one revolution to be 

r+ - t_ = 2.066 x 10~ 10 sec (77) 

This is a small effect that vanishes when the origin of Fermi coordinates moves along a straight line or when w = 0. 



VII. COORDINATE SPEED OF LIGHT 



The three dimensional space of Minkowski space-time is homogeneous and isotropic, and consequently, the speed of 
light is a constant. However, for the case of Fermi coordinates of an observer moving in a circle, the symmetry of the 
3-d space is broken by an acceleration. Consequently, the space is not homogeneous and isotropic and the coordinate 
speed of light depends on position in the 3-d space. 

For a general space-time metric gtj , the coordinate speed of light v c in the direction of a unit vector e Q is given by 
@ 



where 



Vc = TT a C ( 78 ) 
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and the unit vector satisfies the relation e a e@ = 1, where the 3-d spatial metric is given by 42 



90a gop , orA 

7a/3 = 9af3 ~ (80) 

Equation ( |78| ) shows that for a general space-time metric the speed of light depends on direction e a and on position 
and time, through the metric components. A non-zero off-diagonal metric component go a leads to an anisotropy 
(directional dependence) of the coordinate speed of light. The metric for Fermi coordinates given in Eq. ( p4[ ) has 
vanishing off-diagonal terms, 7 Q — 0, and consequently the coordinate speed of light is isotropic. However, the metric 
component G( o) has a non-trivial position and time dependence, and Eq. (|7^) gives the coordinate speed of light: 

v c = c(l + (81) 

where C is a function of Fermi coordinate position and time, given by Eq. (|56|). Therefore, I can write the coordinate 
dependence of the speed of light as 

— = -- \XW cos(-A(°>) + sm(-X<°))l + 0(z/ 4 ) (82) 
c a L c c - 

This variable speed of light is a kinematic effect due to the non-inertial nature of the Fermi coordinate system. Note 
that at any given position X = {X^\ X^ 2 \ X^} the speed of light depends on time in a periodic way, with the 
orbital period of the reference frame, oj. This expression can be rewritten as 



— = a - X-n + 0(i/ 4 ) (83) 
c V c / 

where the unit vector n = (— cos(^A"(°)) , — sin(^X(°)) ,0) is in the observer's orbital plane and points in the 
direction of acceleration, toward the center of the circle. See Figure [| As the Fermi reference frame moves around 
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the circle, the vector n always points toward the circle's center in Fermi coordinates. At any given time X^°\ Eq. 
( |S3| ) shows that to second order in v, in the plane given by 

X ■ n = (84) 

the speed of light does not differ from it's vacuum inertial frame value c. In the 3-d half-space containing positions 
X closer to the center of the circle, the speed of light is increased. For the 3-d half-space containing positions X that 
are further from the center of the circle, the speed of light is decreased. As time X^ increases, the dividing plane, 
given by Eq. (pi)) , rotates at angular frequency u> around the circle, with the origin of Fermi coordinates remaining 
the point at which the plane is tangent to the circle. 

The magnitude of the inhomogeneous variation in the speed of light in Eq. (|||) is small. For the kinematics of 
GPS, using the numbers given in Table 1, with the coordinates X^ x \ X^ ~ b, so the effect has a magnitude 

v 2 

— b = 1.755 x 1CT 12 (85) 
a 

and has a complicated position and time dependence given by Eq. (H|). Recently, GPS data have been used to obtain 
an upper bound on the anisotropy of the speed of light to be Ac/c < 10~ 9 , see Ref. J43| and references contained 
therein. Given the present accuracy of GPS, the small inhomogeneity of the speed of light, given by Eq. ([S2|)is 
probably not measurable. 



VIII. SUMMARY 



I have considered the problem of an observer moving in a circle around a central point in flat space (no gravitational 
fields). The observer carries with him three gyroscopes and uses them to define a set of non- rotating axes. He keeps 
his axes at a constant orientation with respect to the gyroscopes' spin axes. The observer's axes are Fermi transported 
along his world line and Fermi coordinates are defined. The Fermi time coordinate is just the observer's proper time, 
according to the construction shown in Figure [ij Fermi coordinates play the role of ECI frame coordinates in GPS. I 
investigate the kinematic effects seen on ideal clocks from the point of view of these Fermi coordinates. Specifically, 
I consider two cases: stationary clocks (with respect to Fermi coordinates) and satellite clocks in circular orbit with 
respect to Fermi coordinates. The proper time interval between two events on the stationary clock, given by Eq. 
(p3|), is a complicated periodic function of the clock's Fermi coordinate time and coordinate position in the plane 
perpendicular to the rotation axis. I considered next a clock that orbits the Fermi coordinate origin (similar to a 
satellite) in the equatorial plane of rotation and computed the relation of proper time to coordinate time. I find that 
the proper time on the clock has time dilation (constant rate offset) plus a periodic dependence on Fermi coordinate 
time at the difference frequency, f2 — u>, between satellite clock orbital frequency, Q, and Fermi coordinate frame 
orbital frequency, lu . 

Next, I looked at the time difference of two counter-orbiting clocks in the equatorial plane of rotation, one west- 
moving and the other east-moving. After one revolution, the difference in proper time on the clocks shows a small 
Sagnac-like effect, increasing with the angular frequency of revolution of the Fermi frame, u, given by Eq. ( |76| ) . 

Finally, I computed the coordinate speed of light in the Fermi coordinates. The speed of light is isotropic (same 
in all directions), but is not homogeneous. Specifically, the speed of light depends on Fermi coordinate position and 
time, but not on direction. Equation (83) gives the change in the speed of light from c, its vacuum value in an inertial 
frame. 

I have computed the magnitudes of these effects for numbers relevant to GPS satellite orbits and orbital frequency. 
The effects are all small for GPS, except for the 12 ns periodicity for an orbiting clock, given in Eq. (|67|), that occurs 
at the difference frequency, Q — u>, between satellite orbital frequency £1 and Earth orbital frequency u>. This effect 
is reminiscent of the periodic effect in GPS described in the introduction. However, the effects here are all kinematic 
since gravitational fields are neglected. The inclusion of gravitational fields is expected to significantly modify the 
results. In particular, from the work of Ashby and Bertotti, this 12 ns periodic effect is expected to cancel when 
gravitational fields are included P]28]] . Work is in progress to extend these results to include the effects of gravitational 
fields. 

While the effects discussed in this work are small for the present GPS, these effects are quite significant for future 
navigation and time transfer systems that are expected to have sub-nanosecond time transfer accuracy 
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APPENDIX A: CONVENTIONS AND NOTATION 



I use the notation and conventions of Synge |F2] ], except that my space-time indices on coordinates, x l , take values 
i = 0,1,2,3. Specifically, I use the convention that Roman indices take the values i = 0,1,2,3, and Greek indices 
a = 1, 2, 3. I use the summation convention for repeated indices, when the same index appears in a lower and upper 
position then summation is implied over the range of the index. 

If x l and x l + dx l are two events along the world line of an ideal clock, then the square of the proper time interval 
between these events is dr = ds/ c, where ds is given in terms of the space-time metric as 

ds 2 = -gij dx l dx J (Al) 

and I choose to have the signature +2. When is diagonalized at any given space-time point, the elements can 
take the form of the flat space Minkowski metric given by goo = — 1, g a p = 5 a p, and go a = 0. 

In this work, I take the global inertial frame coordinates to have the Minkowski metric. The observer's frame is 
described by Fermi coordinates, X^ a \ a = 0, 1, 2, 3, with line element 

ds 2 = -G {ij) dX (i) dX^ (A2) 
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FIG. 1. The observer's circular 3-d path is shown against the background of the global inertial Minkowski coordinates x % . 
The position of the observer (origin of Fermi coordinates) and the orientation of the Fermi coordinate axes are shown 

schematically at a given time. 
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FIG. 3. The global inertial coordinate axes x l and Fermi coordinates X' a ' are shown in relation to the coordinate position 
X, unit vector n, and the plane X ■ n = 0, at a given coordinate time X^°\ 



TABLE I. Numerical values of constants. 



constant definition value 

a Earth orbital semi-major axis 1.496X10 11 m 

uo Earth orbital angular velocity 1.99238 xl0~ 7 sec -1 

b GPS satellite semi-major axis 2.656177xl0 7 m 

n GPS satellite angular velocity 1.45842 x 1CT 4 sec" 1 

c vacuum speed of light 2.997924X 10 8 m/s 
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